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ABSTRACT 

Using the standard equation for the slowdown of a neutron star, we derive a formula 
for the braking index via integration rather than the conventional differentiation. The 
new formula negates the need to measure the second time derivative of the rotation 
frequency, v. We show that the method gives similar braking indices for PSR B1509- 
58 and the Crab pulsar to those already in the literature. We point out that our 
method is useful for obtaining the braking indices of moderate aged pulsars without 
the need for long, phase-connected timing solutions. We applied the method to 20 
pulsars and discuss the implications of the results. We find that virtually all the 
derived braking indices are dominated by the effects of (unseen) glitches, the recovery 
from which corrupts the value of v. However, any real, large, positive braking index 
has implications for magnetic field decay and offers support to recent models of pulsar 
evolution. 
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1 INTRODUCTION 

Neutron stars are powered by rotational kinetic energy and 
lose energy by accelerating particle winds and by emitting 
electromagnetic radiation at their rotation frequency, v. The 
rotation frequency thus decreases with time and this slow- 
down is usually described by the relation 



-Kv 



(1) 



Here, if is a positive constant which depends on the moment 
of inertia and the magnetic dipole moment of the neutron 
star and n is the braking index. Conventionally, the braking 



index it derived by differentiation of equation 1, yielding 



(2) 



In a highly simplified model in which the spin-down torque 
arises from dipole radiation at the rotation frequency, one 
expects n = 3. 

Only 4 pulsars have had their braking indices measured 
and all have n < 3. The Crab pulsar has a value 2. 509 ±0.001 



(Lyne, Pritchard & Smith 19881; 


Lyne, Pritchard & Smith 


1993;), PSR 1509-58 has a braking index 2.837±0.001 ([Kaspi 
et al. 1994, PSR 0540-69 has n = 2.04 ± 0.02 (Manch- 


ester & Peterson 198£ ; 


Nagase et al. 199C; Gouiffes, Finley & 



Ogelman 199£ ) and, finally , the Vela pulsar has the low value 



1.4 ± 0.2 ( [Lyne et al. 1996[ ). Melatos (1997) has shown that 
a modification of the simple model, which involves treating 
the neutron star and the inner magnetosphere as one entity, 



allows him to derive values of braking index very close to 
those observed (except for the Vela pulsar). 

Braking indices are very difficult to measure in all the 
other pulsars. For a typical 'old' pulsar with v — 1 Hz, 
i>= 10 -15 Hz/s, the expected v from equation 2 is only 
~ 10~ 30 Hz/s 2 , much too small to measure even over hun- 
dreds of years of timing - the second derivative only con- 
tributes one extra phase rotation every 600 yr! There are a 
number of pulsars with ages ~20 kyr for which one might 
expect to be able to measure n. However, two different ef- 
fects tend to dominate the value of v over that expected from 



spin-down alone. F irst, these young pulsars glitch often ( She- 



mar fc Lyne 1996 ). These glitches lead to discontinuities in 
both v and v making it very difficult to phase connect (i.e. 
count the exact number of rotations of the pulsar) over the 
glitch. Furthermore, the recovery from a glitch can last many 
hundreds of days and the measurement of i> reflects the re- 
covery rather than the intrinsic spin-down. Finally, young 
pulsars have large random variations in arrival times known 
as 'timing noise'. Cordes & Helfand (1980) recognised that 
the timing noise dominates over the intrinsic v by a factor 
~100 in these pulsars and that many early published values 
of braking indices, based on v, were spurious. 



2 BRAKING INDEX BY INTEGRATION 

Instead of differentiating equation 1, we integrate from a 
time t to t + T to obtain 
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l - n 



= —KT = 



from which follows 



l-n 
Hence 



But, (yi/v2) n —i'i I i'i and so 



1 + 



V\ U 2 —V2 V\ 



(3) 



(4) 



(•») 



(6) 



U X V 2 T 

This allows the braking index to be computed without 
the need to measure v. The advantage of this method is that, 
in principle, v and v can be measured over a short interval 
of time and then re-measured 20 yr later without the need 
for a phase connected solution over the whole 20 yr time 
span. 
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Figure 1. Braking index for the Crab pulsar as a function of 
time over the period 1982-1987. The braking index is computed 
using equation ^| above and parameters for v and v from Lyne et 
al. (1988). 2278 data points are shown. 



2.1 Error analysis 

We can rewrite equation li as 



1 \Vl Vl 



(T) 



Thus, for n ~ 3, the value of the expression inside the brack- 
ets must be ~ 2T. Typically, the fractional error in v is much 
larger than that in v. Let h\ and e 2 be the error in V\ and 
V2- The absolute error, E, on the braking index is thus given 
by 



E ■ 



T v 



+ e 2 



(8) 



Taking a young pulsar such as the Crab as an example, 
the parameters for which are v ~ 30 Hz, v ~ — 3.8 x 
10~ 10 Hz/s and e ~ 10~ 15 Hz/s over one month ( [Lyne, 
Pritchard & Smith 1988). If measurements are made ~1 yr 
apart then the error in the braking index is ~0.01. 

For older pulsars, those near 100 kyr, one can obtain a 
reasonably low error on the braking index, provided the time 
span is large enough. For example, PSR B0540+23 was dis- 
covered more than 25 yr ago. Its parameters are v ~ 4 Hz, 
13 Hz/s and e 



2.5 x 10" 



10" 



Hz/s (Arzou 



manian et al. 1994). Equation ^| yields an estimated error in 
the braking index of ~0.1, accurate enough to distinguish 
between various slow-down models. 



We can do this for every pair of values in the table (i.e. 2278 
pairs in total). The results are displayed in Fig. 1. 

It can clearly be seen that for the shortest intervals, 
both measurement error and short-term timing noise dom- 
inate the value of the braking index. For time intervals 
greater than ~500 days, the braking index is stable near 2.5. 
This is consistent with the ~600 days periodicity seen in the 
timing noise by Lyne et al. (1988). The weighted average of 
all the measurements is 2.502 (however, the braking indices 
are not truly independent); the braking index measurement 
from the largest time interval is 2.516 ± 0.003. Both these 
measurements are consistent with the Lyne et al. (1988) re- 
sult. It is not possible to compute a braking index for the 
Crab over a longer period of time. Lyne et al. (1993) have 
shown that the glitches permanently alter the value of v; 
however the inter-glitch value of the braking index remains 
roughly constant. 

Kaspi et al. (1994) published a braking index of 2.837± 
0.001 for PSR B1509-58 based on a phase-connected solu- 
tion across 11 yr of timing data. This pulsar has not been 
observed to glitch in that interval, and the timing noise is 
surprisingly low for such a young pulsar. We applied our 
method by finding a local fit to v and v over a period of 3 
months in 1993 and 1997. Over this interval of 1700 days we 
find n = 2.80 ± 0.03 consistent with the Kaspi et al. (1994) 
result but without the need for a phase connected solution. 



3 DISCUSSION 

3.1 Pulsars with known braking index 

For the Crab pulsar, Lyne et al. (1988) have tabulated u, v 
and their associated errors at monthly intervals over a period 
of five years. In order to determine the braking index, they 
first compute v by performing a straight line fit to the v 
versus time graph. They derive a braking index of 2.509 ± 
0.001 based on the value of v over this interval. We can 
use their tabulated values and derive a braking index from 
equation @ above, without the need to compute v explicitly. 



3.2 Young, Vela-like pulsars 

There are 21 pulsars with | v \ > 10~ 12 Hz/s. Of these, 
4 have their braking indices measured (see above), and 11 
of the remaining 17 have had one or more glitches and are 
expected to glitch every ~10 yr or so (see e.g. Shemar & 
Lyne 1996). For the Vela pulsar, for example, the intrinsic v 
is completely dominated by the recovery of v following the 
regular glitches. Lyne et al. (1996) solved this problem by 
taking v at a fixed interval after each glitch and deriving v 
and hence the braking index from the changing v. Presum- 
ably, if such a method is valid, one could apply equation 
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over the largest time interval available (25 yr) to obtain an 
error on the braking index of ~0.1. 

The six pulsars which are not known to have glitched are 
PSRs B1951+32, B1853+01, B1930+22, B1643-43, B0906- 
49 and J0631+1036. Taking PSR B1951+32 as an example, 
Foster et al. (1994) have shown that this pulsar suffers from 
excessive timing noise. By fitting for higher orders in v they 
obtain a value of v of 2.7 x 10 -22 Hz/s 2 . This is about a 
factor ~50 higher than that expected from the simple spin- 
down, making it virtually impossible to compute a braking 
index. PSR B0906-49 was discovered 12 yr ago, and has 
not glitched since discovery. Using recent timing data, the 
value of v from timing noise is ~ 1.8 x 10 -23 Hz/s 2 , more 
than 30 times the expected spin-down value. Thu s, timing 
noise is dom inating the value of v in the literature (D'Amico 
et al. 1988) and a braking index cannot be computed. In 
particular the formal fitting errors given by TEMPO (the 
least-squares pulsar timing package) are seriously underes- 
timating the systematic effects of the timing noise. 

For these young pulsars, the new method suffers from 
similar problems to the standard calculation of the braking 
index. In Vela-like pulsars glitches dominate the overall spin- 
down behaviour, masking the true spin-down v. In slightly 
older pulsars, glitches are less frequent; however, the true v 
is hidden in large timing noise residuals. 



3.3 Pulsars with moderate v 

As described above, it should be possible in principle to 
obtain an error on the braking index as low as ~0.2 for 
pulsars of age ~ 10 5 yr discovered more than 20 yr ago 
(especially those with low timing noise). We select pulsars 
in the following way. We construct an 'expected' value of v 
for every pulsar, assuming a braking index of 3 and select all 
those with v> 10~ 27 Hz/s 2 . We reject all those pulsars which 
are known to have glitched in the past and the 21 pulsars 
described in 3.2 above. We then searched the literature for 
at least two timing solutions for the remaining pulsars and 
with a v of sufficient accuracy to allow the error on the 
computed braking index to be less than 20. We discovered, 
however, that in some publications, the epoch of v and v 
corresponded not to the middle of the data span, but rather 
to the start. If there is a significant value of v in the data 
span which is unaccounted for (either related to timing noise 
or to the intrinsic spin-down), this leads to an error in the 
quoted v of ~ | i> T which, for large T, can be significantly 
larger than the formal fitting error. This led to the rejection 
of a further five pulsars. 

Table 1 gives the braking index for 20 pulsars which 
survive the above selection criteria. Column 1 of the table 
gives the pulsar name. Columns 2 and 3 give the pulsar's 
rotation frequency and its first derivative, and column 4 the 
time interval between the epochs of measurement. Column 
5 gives the braking index and associated error according to 
equations ^ and |§| The final column gives the references 
for the timing solutions. Note that for PSR B0656+14 the 
value of v given in Ashworth & Lyne (1981) is clearly in 
error. Hence, the time interval between timing solutions is 
rather shorter than it might have been. 

The values of the braking indices listed in the table 
raise the question as to how realistic the quoted errors are. 
Typically (but not always) in timing solutions the errors are 



twice the formal standard error given by the TEMPO least- 
squares timing package. However, the errors are computed 
assuming that the each TOA has uncorrelated residuals (i.e. 
white noise), which is clearly not the case in pulsars with 
significant (red) timing noise. Thus the error on v is likely 
to be larger than that given from the fit alone. Also, as 
described above, not fixing the epoch in the middle of the 
data span leads to errors in the value of v. 





PSR B 


V 

(Hz) 


v xl0~ 15 
(Hz/s) 


Interval 
fdavs) 


71 


Refs 


1 


0114+58 


9.86 


-562.3 


2271.1 


-9.6+1.5 


1,2 


2 


0136+57 


3.67 


-144.2 


4492.0 


-81+4.7 


5,2 


3 


0154+61 


0.43 


-34.2 


4336.5 


28+14 


5,2 


4 


0540+23 


4.06 


-255.0 


5543.5 


11.1+8.6 


8,2 










5990.5 


11.81+0.12 


3,2 


5 


0611+22 


2.99 


-530.8 


5541.5 


20.1+1.1 


3,2 


6 


0656+14 


2.60 


-371.5 


2163.3 


14.7+1.4 


1,2 


7 


0740-28 


6.00 


-604.9 


4245.2 


17.7+1.4 


9,7 










5827.2 


25.6+0.8 


9,2 


8 


0919+06 


2.32 


-74.0 


4521.7 


28.9+4.1 


5,2 


9 


1221-63 


4.62 


-105.7 


6661.3 


18.7+12.3 


6,10 


10 


1356-60 


7.84 


-389.9 


7050.3 


6.3+6.8 


6,10 


11 


1719-37 


4.23 


-194.6 


4824.0 


-183+10 


6,2 


12 


1742-30 


2.72 


-79.0 


1581.0 


-132+5 


7,2 


13 


1829-08 


1.54 


-151.4 


1541.0 


2.5+0.9 


4,2 


14 


1907+10 


3.52 


-32.8 


5842.5 


24+17 


3,2 


15 


1914+09 


3.70 


-34.5 


5556.6 


-15+16 


3,2 


16 


1915+13 


5.14 


-190.0 


6080.5 


36.08+0.48 


3,2 


17 


2000+32 


1.44 


-216.7 


1381.0 


-226+4.5 


4,2 


18 


2002+31 


0.47 


-16.7 


6076.5 


23.3+1.0 


3,2 


19 


2148+52 


3.01 


-91.2 


2307.2 


49.6+3.5 


1,2 


20 


2334+61 


2.02 


-776.2 


2347.1 


8.60+0.13 


1,2 



[Table 1 references: 1-Dewey et al. (1988), 2-Arzoumanian 
et al. (1994), 3-Gullahorn & Rankin (1978), 4-Clifton et al. 
(1992), 5-Backus, Taylor & Damashek (1982), 6-Newton, 
Manchester & Cooke (1981), 7-Siegman, Manchester & Dur- 
din (1993), 8-Helfand et al. (1980), 9-Manchester et al. 
(1983), 10-Parkes timing programme.] 

Figure 2 displays the braking index for these 20 pul- 
sars with error bars a factor of 5 larger than in the table. 
We believe these error bars are conservative and more accu- 
rately reflect the contributions from the underlying timing 
noise to the intrinsic value of v. From the figure, 14 of the 20 
pulsars have significant values of braking index. The eight 
pulsars with the smallest error bars all have positive brak- 
ing indices (PSRs B0540+23, B0611+22, B0656+14, B0740- 
28, B1915+13, B2002+31, B2148+52 and B2334+61). The 
four pulsars with large negative braking indices are PSRs 
B0136+57, B1719-37, B1742-30, B2000+32 and B2255+58. 



3.4 Implications 

The presence of glitches in pulsars can lead to spurious 
values of the braking index and we surmise that (unseen) 
glitches are the main cause for n ^ 3. If we assume that 
glitches cause A v / v to be positive in all cases and that 
the value of v recovers to nearly its original value, then large 
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Pulsar Number 



Figure 2. Braking index for 20 pulsars from Table 1. The error 
bars are a factor of 5 larger than in the table. 



negative braking indices can be caused by glitches between 
timing solutions and positive braking indices arise when the 
glitch occurs before the epoch of the first timing solution. Al- 
par & Baykal (1994), in a statistical survey of the frequency 
of pulsar glitches, showed that the number of glitches, g, in 
a sample of pulsars is given by 



(9) 



where T is the time interval ov er which the pulsar is moni- 
tored and 8u/u is 1.74 x 10" 4 (lAlpar & Baykal 19941). 



Using the values of v, v and T from the table, we derive 
g = 2.8, i.e. one expects there to have been glitches in ~3 
pulsars in between their timing solutions. We have four pul- 
sars with large negative braking indices which is consistent 
with this estimate. The values of A v / v for these four pul- 
sars are 1.3, 3.6, 0. 5, 4 and 0.6 x 10~ 3 r espectively, values 
typical of glitches ( Shemar fe Lyne 1996] ). It is thus possi- 
ble that these pulsars have all glitched between the epochs 
of the observations. However, none of the 4 appear to have 
suffered any change in v, i.e. any step change at the time 
of the glitch would have been recovered in the intervening 
time period. Alternative explanations are hard to find - for 
PSR B2000+32, for example, one of the two timing solu- 
tions would have to be in error by more than 100a to obtain 
a 'real' braking index of 3. 

For the positive values of braking index, it is harder to 
compute the number of glitches expected before the timing 
solutions. Both the recovery time and the form of the recov- 
ery (linear or exponential) for the change in i> are unclear. 
Shemar & Lyne (1996) estimate that the recovery time is 
greater than 3 yr for the exponential case or tens of years if 
the recovery is linear. Indeed, in some pulsars, a permanent 
alteration in the value of i> is observed. If we let T = 25 yr, 
then we would expect 6.4 glitches prior to the first timing 
solution for these 20 pulsars. We have 8 positive braking in- 
dices, broadly in line with this estimate (and we also note 
that the 2 pulsars with the highest glitch probability are 
PSRs B2334+61 and B0611+22, both of which have posi- 
tive braking indices). 

A number of other factors may contribute to n / 3. 
It is possible that the different methods of obtaining the 



pulsar position could affect the results (however, for PSR 
B0540+23, for example, the positions are identical within 
the errors). Another possibility is that the different fitting 
software, clock corrections, ephemeris changes etc can affect 
the result but it is difficult to believe this could have an effect 
at the 500cr level. These possibilities could be overcome by 
recovering the original TOAs from the 1970s and reducing 
them with the same software as the 1990s data. 

The glitch interpretation put forward above, is statisti- 
cal in nature. It may be that in one or more cases the braking 
index is indeed significantly larger than 3. What are the im- 
plications of this? Blandford & Romani (1988) showed that 
if the value of K in equation 1 is time variable then one can 
write 



K v 
K v 



(10) 



K/ K then gives the timescale for variation, if we assume 
that n = 3. For the cases in which n b s > 3, K< 0, and the 
timescale is surprisingly constrained in the range ~14 to ~56 
kyr. K< can arise from magnetic field decay, alignment of 
the spin and magnetic axis, or an increase in the moment of 
inertia. Evidence for and against magnetic field decay and 
axis alignment have raged in the literature for 25 yr without 
any clear consensus emerging. Recently, Chen, Ruderman & 
Zhu (1998) have developed a model in which the magnetic 
field increases by a factor of ~10 over the first 10 kyr of a 
pulsar's life (which may explain the low braking index of the 
Vela pulsar) and then decreases thereafter. In their model, 
they expect a braking index of >5 for middle aged pulsars 
before it settles in older pulsars back to 3. If any of the above 
braking indices are thus real, this implies a magnetic field 
decaying by a factor of ~100 in ~200 kyr, lending support 
to the Chen et al. (1998) model. 

In any case, one is forced to conclude that the initial 
starting premise that the spin-down of a pulsar can be de- 
scribed as a simple power law in v is highly questionable. In 
virtually all young pulsars, and in most older ones, complex 
glitch behaviour appears to be dominating the spin-down 
behaviour. This glitch behaviour, coupled with changes in 
the magnetic field strength over time, indicates that the en- 
tire concept of a 'braking index' must be treated with some 
caution. 

Verification of the glitch hypothesis or the presence of a 
true braking index could be obtained by deriving additional 
timing solutions for these pulsars now. This would give a 7 
yr span back to the Arzoumanian et al. (1994) data and, in 
some cases, more than 20 yr back to earlier data. One might 
then expect to observe the exponential recovery of v in the 
case of a glitch. 



4 CONCLUSIONS 

We have derived a new method for computing pulsar brak- 
ing indices based on integration of the standard slow-down 
equation. We show that the method replicates the known 
values for both the Crab pulsar and PSR B1509-58. The new 
method conveys no advantage over the previous method for 
pulsars with regular glitches and/or large timing noise. The 
advantage of equation ^| is that for pulsars with ages ~ 10 5 
yr where both the timing noise and the glitch frequency are 
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relatively low, it is not necessary to obtain a phase connected 
solution across 25 yr or more in order to obtain a braking 
index. Rather, local fits over 1-2 yr at large intervals may 
be sufficient to yield the braking index to an error of ~0.1. 

We computed braking indices for 20 pulsars based on 
timing solutions available in the literature. Four pulsars have 
large negative braking indices which we believe have been 
caused by glitches occurring in between the epochs of the 
timing solutions. A further eight pulsars have moderate posi- 
tive braking indices which appear to be robust to the effects 
of timing noise but are possibly due to glitches which oc- 
curred before the start of the timing observations. If any of 
the large, positive braking indices are real, it may provide 
evidence for magnetic field decay in moderate aged pulsars. 
However, the entire concept of a smooth pulsar spin-down 
and of a constant braking index must be treated with some 
caution. 
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